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A fundamental relation between the metallicity, gas 
content, and stellar mass of local galaxies 
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ABSTRACT 

Recent results have suggested that the well known mass-metallicity relation has a 
strong dependence on the star formation rate, to the extent that a three dimensional 
'fundamental metallicity relation' exists which links the three parameters with mini- 
mal scatter. In this work, we use a sample of 4253 local galaxies observed in atomic 
hydrogen from the ALFALFA survey to demonstrate, for the first time, that a similar 
fundamental relation (the Hi-FMR) also exists between stellar mass, gas-phase metal- 
licity, and Hi mass. This latter relation is likely more fundamental, driving the relation 
between metallicity, SFR and mass. At intermediate masses, the behaviour of the gas 
fundamental metallicity relation is very similar to that expressed via the star forma- 
tion rate. However, we find that the dependence of metallicity on Hi content persists 
to the highest stellar masses, in contrast to the 'saturation' of metallicity with SFR. 
R is interesting to note that the dispersion of the relation is very low at intermediate 
stellar masses (9 < log (M*/M0) < 11), suggesting that in this range galaxies evolve 
smoothy, in an equilibrium between gas inflow, outflow and star formation. At high 
and low stellar masses, the scatter of the relation is significantly higher, suggesting 
that merging events and/or stochastic accretion and star formation may drive galaxies 
outside the relation. We also assemble a sample of galaxies observed in CO. However, 
due to a small sample size, strong selection bias, and the influence of a metallicity- 
dependent CO/H2 conversion factor, the data are insufficient to test any influence of 
molecular gas on metallicity. 

Key words: galaxies: evolution - galaxies: formation - galaxies: abundances - galax- 
ies: statistics 



1 INTRODUCTION 

Understanding the physical processes regulating the metal- 
licity of galaxies is one of the most important problems in 
modern astrophysics. The gas-phase metallicity functions 
as a fossil record of a galaxy's history, encapsulating in- 
formation about both the past star formation and the in- 
fall/outflow of interstellar gas. 

It has long been known that metallicity correlates 
well with both bolometric luminosity and stellar mass, 
with the latter (the 'mass-metallicity' relation) in partic- 
ular bein g extensively studied in recent years, b oth lo- 
cally (e.g. iLequeux et al.lll979l ; iTremonti et al.ll200 j ) and at 
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higher redshift s (e.g. ISavaglio et"alll2005l: lErb et al.l 120061 



Maiolino et~aT] |200S| ; llvlannucci et all 120091 : ICresci etlil] 
20121 ). A picture has emerged of a well-constrained local re- 
lation which evolves strongly out to z > 3, placing valuable 
constraints on galaxy evolution models. 

A host of physical models has been proposed to ex- 
plain the mass metal licity relation , inclu ding star forma- 
tion e fficiency issues (IBrooks et al.l 120071 ') and IMF varia- 
tions (jKoppen et al.ll2007r ). Perhaps unsurprisingly, though, 
most explanations for the relation have focussed on ex- 
planations involving th e global behaviour of gas (i.e. 
lAudouze fc Tinslevlll976r). which can serve both as a pris- 



tine i nfalling diluent (e.g. iDave et al.ll2010);lFinlator 
2008) and as a metal rich outflow ( Spitoni et al 



Peeples fc Shankarll201ll ). 
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More recently, it was found that the mass-metallicity 
relation can be seen as a consequence of a more fundamen- 
tal relation betwee n stellar mass, metallicity, and star for- 
mation rate (SFR) (IMannucci et al.ll2010l ; lLara-Lopez et al.l 
l20ld ; lAndrews fc Martinil 120131 '). This 'fundamental metal- 
licity relation' ('FMR' hereafter) was found to exhibit re- 
markably little scatter, and - more i mportantly - was found 
not to evolve out to at least z — 2.5 (|Mannucci et al.1120101 ; 
lLara- Lopez et al.l 120101 ). As metallicity anti-correlates with 
SFR at a constant stellar mass, the apparent evolution of 
the mass metallicity relation can be explained as a natural 
consequence of the fact that samples of galaxies at higher 
redshifts have elevated star formation rates relative to local 
galaxies. If the FMR is seen as a relation in a 3D parameter 
space defined by stellar mass, metallicity, and SFR, sam- 
ples of galaxies at different redshifts simply inhabit different 
regions of this unchanging 'fundamental plane'. 

The strength of the correlation between these three pa- 
rameters is at first glance somewhat surprising, and strongly 
suggests the presence o f highl y efficient regulatory mech- 
anisms. IMannucci et al.l (|2010l ) attribute the existence of 
the FMR to a well-regulated balance between pristine in- 
flows and enriched outflows. In this model, metal-poor gas 
falls into a galactic halo, and acts to both boost star for- 
mation and decrease metallicity by a 'dilution' effect. Like- 
wise, star formation-driven outflows remove metal-rich gas 
at a rate which is strongly dependent on stellar mass. These 
two effects combine, resulting in a tight correlation between 
the three parameters stellar mass, metallicity, and SFR. 
These mech anisms have b een fu rther explored bv lDave et al.l 
(|2012f ) and iDaval et al.l (j2012T >. who each present a set of 
simple equations which link the observable quantities of the 
FMR to the behaviour of infalling and o utflowing gas. The 
difference between these models is that iDave et alj (|2012l ) 
assume an 'equilibrium condition', whereby the inflow rate is 
balanced by the sum of the star formation and outflow rates 
(this is some times referred to as the 'reservoir' or 'bathtu b 
model' - see lBouche et al.ll2010l ; lKrumholz fc Dekelll201lh . 

The common thread underpinning all physical and an- 
alytical models of the FMR is the understanding that the 
driving force behind these tight correlations is the be- 
haviour of gas. In particular, as mentioned above, the anti- 
correlation between metallicity and SFR is thought to be a 
byproduct of a more fundamental relation between metallic- 
ity, gas content, and stellar mass. It has already been noted 
by authors that there is a connectio n between metallicity 
and gas content ; lLequeux et alj (|l979l ) found an inverse cor- 
relation between heavy element abun dance and gas f r action 
for l ocal blue com p act g alaxies. Both lSkillman et al.l (|1996T ) 
and lPeeples et alj (|2008l ) noted t hat local metal rich galax- 
ies tended to be gas-poor, and iBoselli fc Gavazzil (| 20061 ) 
made the same obs ervation for cluster galaxies. Recently, 
Irlughes et alj (|2013T ) systematically analysed the influence 
of Hi gas fractions on the mass-metallicity relation for a 
sample of 260 local galaxies, finding that the gas fraction 
is clos ely related to the scatt er in the mass-metallicity re- 
lation. iRobertson et alj (|2012P ) also found a connection be- 
tween gas-phase metallicity and Hi deficiency for a small 
sample of 6 local cluster galaxies, with Hi deficient galaxies 
being more metal rich by ~ 0.1 dex. 

In this paper, we use a large sample of local galaxies to 
reveal the existence of a fundamental relation between stellar 



mass, metallicity, and atomic gas content. We also explore 
the relation between stellar mass, metallicity, and molecular 
gas content; however, such a relation is much more difficult 
to investigate due to the limited data available, selection 
effects, and issues associated with the CO/H2 conversion 
factor. The paper is structured as follows: in §2 we present 
the samples used in this analysis, and describe the method 
used to determine metallicity. §3 gives our results, and we 
discuss the implications in §4. We present our conclusions 



2 SAMPLE SELECTION 

The primary motivation behind our sample selection was 
assembling a large number of galaxies with the full com- 
plement of data, sufficient for measuring the gas mass, the 
stellar mass, the gas-phase metallicity, and the star forma- 
tion rate in a homogeneous manner across the sample. By 
far, the most time consuming of these parameters to obser- 
vationally measure is the gas content. While surveys such 
as the SDSS now routinely measure the optical properties 
of millions of z ~ galaxies, obtaining hydrogen masses for 
similarly large numbers of galaxies remains more of a chal- 
lenge. While large, modern surveys for atomic hydrogen can 
now collect data for (on the order of) 10 4 galaxies, molecu- 
lar hydrogen remains more elusive yet, with fewer than 10 3 
galaxies having the CO observations required to extrapolate 
a molecular gas mass. 

Below, we discuss the surveys we have used to build our 
samples of local galaxies. 



2.1 HI sample 

To construct our Hi sample, we used the first data release 
of the 'Arecibo Legacy Fast ALFA' survey ( ALFALFA), a 
wide- field sur vey of Hi in the local U niverse (JHavnes et alj 
l201ll ; see also iGiovanelli et alj 120051 ). This work uses data 
from the initial 'a. 40' data release, consisting of 15,855 
ga laxies. The details of the survey selection can be found 
in lHavnes et alj|201ll - in brief, however, the current data 
release of the survey covers 2800 deg 2 over two regions of 
the sky with a beam size of 3.5', and selects galaxies out to 
~ 250 Mpc (z ~ 0.05). For such nearby galaxies, the SDSS 
fibre often only covers a fraction of the galaxy; see Appendix 
A for a discussion of potential aperture biases caused by the 
SDSS fibre. With a large number of sources and good sensi- 
tivity, the ALFALFA a. 40 is an ideal resource for examining 
large scale trends in atomic gas content. 

This database of Hi-observed galaxies was then 
cross-matched with the SDSS spectroscopic survey 
IjAbazaiian et alJl2009T ). To cross-match the two databases, 
we defined a search box around each detected Hi position 
of 10 arc seconds in RA and DEC, and a redshift interval 
Az — 0.001. We discarded any sources that returned either 
zero matches, or more than one unique match. The total 
number of galaxies in the ALFALFA-a.40 survey appearing 
in the SDSS spectroscopic database is 8930. 

Of these galaxies, it was necessary that all lines re- 
quired for metallicity diagnostics were detected. We dis- 
card any galaxies with a Ha S/N lower than 25, which 
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Figure 1. The mass- metallicity relation for the 4253 ALFALFA galaxies in our sample. The grey shaded areas show the area that 
contains 64% (light shaded area) and 90% (dark shaded area) of all SDSS galaxies from the Mannucci et al. (2010) study. The coloured 
lines show the mean trends for galaxies in bins of Ho-SFR (left panel) and Hi mass (right panel). It can be seen that at a constant stellar 
mass, metallicity is a decreasing function of both SFR and Hi mass. 



has the effect of ensuring that we have the full comple- 
ment of lines required for our metallicity diagnostics. As 
discussed further in £|2-3I we use two metallicity diagnostics: 
the [Nll]A6584/Ha ratio, and the well-known 'R23' parame- 
ter, defined as R23 = ([OILU3727 + [OIIIJA4958, 5007)/H/3. 

Of the initial 8930 galaxies, 66% have a Hq S/N greater 
than 25, leavi ng us with 5952 galaxi es. 

Following iMannucci et al.l ( 201(3 ). we also discard galax- 
ies for which our two metallicity diagnostics suggest 
metallicities differing by > 0.25 dex, any galaxies which 
are 'AGN dominated' using the BPT cut presented by 
IKauffmann et alj (J2003T I. 

After these cuts, the final Hi sample consists of 4253 
galaxies. This sample is nearly 20 times larger than the sam- 
ple used by Hughes et al. (2012) for a similar study, enabling 
us to achieve unprecedented statistical significance, and to 
investigate a range of galaxy properties unexplored by pre- 
vious studies. 

For all galaxies appearing in the SDSS, total stellar 
masses (i.e ., corrected for the effect of the fibre aperture) are 
taken from IKauffmann et al.l (|2003l ). We have calculated star 
formation rates within the SDSS fibre, allowing direct com- 
parison with Mannucci et al. (2010). Star formation rates are 
therefore calculated using the luminosity of the extinction 
corrected Ha em ission line, usin g the Ha /SFR c onversion of 
iKennicuttJ (|l99cf ). modified to a lChabrierl (|2003T l IMF using 
a conversion factor of 1.58 (this is the same star formation 
conversion adopted by Mannucci et al. 2010). 



2.2 H 2 sample 

The study of the molecular phase of the ISM is more ob- 
servationally challenging than the atomic phase. Molecular 
hydrogen lacks a permanent electric dipole, necessitating the 
use of a proxy molecule (typically CO), from which the mass 



of H2 can be extrapolated. The observation of CO is still 
tricky, however, and extragalactic CO surveys remain less 
populous than Hi surveys by well over an order of mag- 
nitude. The result, of course, is that samples dealing with 
molecular hydrogen observations are statistically impover- 
ished relative to surveys of Hi. 

To minimise this issue, we assemble a sample of galaxies 
from the literature that have both compatible metallically 
diagnostics (i.e. R23 and [NII]/Ha) and CO data available. 
These additional galaxies have physical parameters (such 
as stellar mass) derived differently to the larger samples of 
SDSS galaxies. Due to the varying selection parameters, we 
take care to separately identify all galaxies which are not 
part of the large surveys. 



2.2.1 COLDGASS sample 

The primary source of H2-obse rved galaxies in th is 
work is the COLDGASS survey (|Saintonge et all l201ll ). 
COLDGASS is the H 2 subset of the GASS survey 
(|Catinella et alj|20ich . a study designed to obtain deep Hi 
data for a sample of relatively massi ve z ~ galaxies. The 
initial data release of COLDGASS (|Saintonge et alj|201ll ) 
consists of 307 galaxies observed in CO, with a detection 
rate of 54%. 

After applying the same selection criteria as applied 
to the ALFALFA sample above (i.e. a robust SDSS coun- 
terpart, well-detected metallicity indicators, lack of 'AGN 
domination', and Ha S/N > 25, along with a CO detec- 
tion), we are left with 72 galaxies. However, a further 33 of 
these galaxies are have highly inconsistent metallicity indi- 
cators; the [Nll]A6584/Ha ratio exceeds the 'R23' parameter 
by > 0.5 dex (in fact, in all of these cases the metallicity de- 
rived from the [Nll]A6584/Ha ratio simply saturates at the 
maximal value allowed by the calibrations). These galax- 
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Figure 2. The mctallicity-SFR (left panel) and the metallicity-Hl (right panel) relation. The coloured lines show the mean trends for 
galaxies in bins of stellar mass. In each plot, a systematic dependence of metallicity on both SFR and Hi mass can be seen. At the 
highest stellar masses there is little dependence of metallicity on SFR, while the dependence of metallicity on Hi mass is seen across the 
full range of the sample. 



ies most likely host LINERs or some weak Seyfert nucleus 
that contaminates the emission lines. In either of these cases 
the metallicity inferred from emission lines is not reliable. 
After excluding the COLDGASS galaxies with saturated 
[Nn]A6584/Ha metallicities, we are left with 39 COLDGASS 
galaxies observed in CO. 

Galaxies in the COLDGASS sample were cross-matched 
against the SDSS database as described above, which again 
provided star formation rates and stellar masses, just as for 
the ALFALFA sample. 



2.2.2 Other CO galaxies 

As explained above, the number of galaxies currently ob- 
served in CO is approximately two orders of magnitude 
smaller than those observed in Hi. As a result, our H2 sam- 
ple is necessarily far smaller than our atomic gas sample. To 
go some way towards rectifying this, for the purposes of this 
study we also include a number of local galaxies that have 
been observed in CO, but are not part of a single homoge- 
neous survey (like COLDGASS). Note that we take care to 
identify these galaxies separately in all plots hereafter, in 
order to avoid confusion. 

We include an additional 23 literature galaxies 
observed in CO. These galaxies were prima ril y as - 
sem bled from the literatur e by ISchruba et al.l (120121) 
and lObreschkow fc Rawlingsl (|2009l ). While ISchruba et al.1 



(|2012r i give metallicities taken from a variety of sources (pri- 
marilv lMoustakas et al.ll2010T ). in order to remain as consis- 
tent as possible we re-derive metallicities using the method 
described above (i.e., taking the mean of the metallicity 
derived using the [Nll]A6584/Ha ratio and R23). Where 
the full complement of lines were unavailable, making this 
method impossible, we adopted the metallicity given by a 



single diagnostic - either R23 or the [Nll]A6584/Ha ratio. 
This was only the case for 2 of these galaxies, however. 

Most of the galaxies asse mbled by Schruba et al.l (|2012T l 
appear in the SINGS survey (|Kennicutt et al.ll2003l ). They 
have stellar mas ses derived using SED fitting, given by 
iNoll et al.l (|2009l f1. and Hi masses from T he Hi Nearby 
Galaxy Survey, THINGS (jWalter et al.ll2008l ). 

5 of our supplemental CO galaxies did not have SED- 
derived stellar masses available - stellar masses for these 5 
galaxies were derived using K-hanA magnitudes (as listed 
in the NASA ADS), convolved with a {B — V) colour- 
appropri ate mass-to- light ratio , following the stellar models 
given bv lBell fc de Jond (|200ll ). In addition, 7 of our supple- 
mental CO galaxies did not have a THINGS Hi mass listed 
- for these galaxies, Hi masses were derived using the 21 cm 
magnitude taken from the online database HYPERLEDA. 

To summarise, the local molecular gas sample consists 
of 62 galaxies, 39 of which are taken from the homoge- 
nous, SDSS-crossmatched COLDGASS database, with phys- 
ical parameters derived identically as for the ALFALFA-Hl 
sample described above. To improve our sample size, we 
have added 23 local galaxies taken from smaller CO surveys. 
These secondary galaxies have physical parameters derived 
using different methods to the SDSS samples, and therefore 
do not constitute a homogeneous dataset. 

The advantage of CO observations, in contrast with ob- 
servations of atomic hydrogen, is that it can be detected out 
to cosmological distances with current instrumentation. The 
detection of Hi at redshifts beyond 2 ~ 0.3 will not be pos- 
sible until the commissioning of the Square Kilometer Array 
(see i|4]below). In order to examine the relationship between 
H2 content and metallicity at high redshift, we also include 
a sample of high-2 star forming galaxies which have well 



1 We have scaled the given stellar masses to a Chabricr IMF. 
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Figure 3. The mean scatter in the 'fundamental metallicity relation' as a function of the (SFR) parameter a (left panel, and Eq. 2) 
and the (Hi) parameter /3 (right panel, and Eq. 3). Solid lines show the dispersion defined as the scatter around the polynomial best fit 
to the (a = 0) M*-Z relation, taken from Mannucci et al. (2010). Dashed lines show the dispersion defined as the scatter around a best 
fitting fourth-order polynomial, recalculated at each a. In each case, the dispersion in the mass-metallicity relation can be reduced by 
including an additional dependence on either SFR or Hi mass. The minimum dispersion reached for the SFR and Hi relations is 0.073 
dex and 0.075 dex (respectively). We estimate that the uncertainly on the metallicity is ~ 0.05 dex (excluding systematics). 



characterised stellar masses, CO masses, and metallicities. 
We use a sample of 'sub-millimeter galaxies' (SMGs), which 
have metallicities est imated using their [N u]A6584/Ha ra- 
tio (using data from Swinba nk et al.l 120061 ) . molecular gas 
masses tak en fromiBothwell et alj ( 20131 ) . and stellar masses 
derived bv lHainline et al.l ( 20111 ). 



2.2.3 H2 gas masses 

The CO-to-H2 conversion factor {aco = Mh2/L^ O ( 1 _ j), is 
one of the main sources of uncertainty in determining the 
mass of molecular gas from the CO luminosity. The main 
issues involved in the determination of aco , and its depen- 
dence of the various galaxy parameters, are reviewed in Bo- 
latto et al. (2013). For our study, a central concern is the 
strong dependence of the conversion factor on metallicity, 
which can obscure trends between H2 mass and metallicity. 
We discuss this issue further in H3.3I below. In this work, 
we have adopt a metallicity-dependent conversion factor, as 
given by Schruba et al. (2012): 



ratios involving bright lines such as [OIII]5007, [Oil] 3727, 
[NII]6584 and the Balmer lines, are sensitive to the metal- 
licity, although, in many cases with large scatter. The cali- 
bration of these 'strong line ratios' has been attempted both 
using empirical methods (e.g. by exploiting 'direct' metallic- 
ity tracers) and through photoionisation models, or through 
mixed methods (see discussion in Kewley & Ellison 2008 
and in Maiolino et al. 2008) . Each of these methods has ad- 
vantages and disadvantages, but it is beyond the scope of 
this paper to discuss the various metallicity calibrations in 
detail. Since in this paper we wish to compare our results 
with the FMR inferred by Mannucci et al. (2010), here we 
adopt the same calibration adopted by those authors: i.e. 
the strong line calibrations given in Maiolino et al. (2008), 
which use a combination of models and empirical calibra- 
tions. In particular, we use the R23 parameter and [Nil] /Ha 
ratio as metallicity diagnostics, and we follow the same cri- 
teria as Mannucci et al. (2010), in that we take as the final 
metallicity the mean of the metallicities determined using 
the two methods. 



log aco = log(A) + JV(12 + log O/H) - 8.7, 



(1) 



with A = 8.2 ± 1.0 and N = -2.8 ± 0.2. This form for 
aco is very similar to that given in other recent works (e.g. 
iNaravanan et al.ll201ll ; lGenzel et al.ll201ll ), especially at the 
relatively high metallicities (12 + log O/H > 8.5) probed by 
our sample. 



2.3 Metallicity determination 

The determination of the gas metallicity through strong line 
ratios has been widely discussed in the literature. Some line 



3 RESULTS 

3.1 The fundamental metallicity relation and HI 
content 

Figure [T] shows the mass metallicity relation for our sample 
of 4253 local galaxies appearing in both the ALFALFA a. 40 
catalogue, and the SDSS. The grey shaded areas shows the 
region containing 64% and 90% (light grey a nd dark grey re- 
specti vely) of all SDSS galaxies, as given bv lTremonti et al.l 
(|2004h . 

The coloured lines in each panel show the mean metal- 
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licity of ALFALFA galaxies, separated into bins of SFR (left 
panel) and Hi mass (right hand panel). Bin sizes are 0.1 dex 
in Hi and 0.15 dex in SFR. Bin sizes were chosen to en- 
sure sufficient numbers of galaxies in each bin; the number 
of galaxies per bin ranges from 21 (at the high end of SFR 
and Hi mass) to 841. 

The left hand panel, displaying the mass-metallicity re- 
lation separated into bins of SFR, shows data which is func- 
tionally equivalent to the result presented by Mannucci et al. 
(2010) - albeit with a smaller sample size. The global trend 
reported by those authors is clearly evident, with metallicity 
increasing with stellar mass, and decreasing with SFR at a 
constant stellar mass. 

In the right panel of Fig. [T] it can immediately be seen 
that there is also a strong, systematic influence of Hi mass 
on the mass metallicity relation, in much the same way as 
emerges with the SFR. At a given stellar mass, galaxies with 
greater Hi masses have proportionally lower metallicities. 
This effect can be seen down to the lowest stellar masses 
probed (log(M*/Mo) ~ 8.7) - below which, the number 
of galaxies per bin becomes small, and any trends become 
unreliable) . 

Interestingly, the two relations show markedly different 
behaviour at the highest stellar masses. As shown in Fig. 
1, and as noted by Mannucci et al. (2010), the correlation 
between metallicity and SFR weakens towards the highest 
stellar masses, and the most massive galaxies in the sample 
(log (M,/Mq) > 10.5) show no correlation between metal- 
licity and SFR. Conversely, the correlation between metal- 
licity and Hi mass (as can be seen in the separation of the 
various Hi mass bins in metallicity space) is evident across 
the full range of stellar masses. 

This can be better seen in Fig. [2] which shows metallic- 
ities plotted as a function of both SFR (left panel) and Hi 
mass (right panel). The mean metallicity is shown for bins 
of varying stellar mass. Less massive galaxies show a strong 
inverse dependence of metallicity on SFR, while the most 
massive galaxies - log(M,/MQ) > 10.5 - show an approx- 
imately constant metallicity with SFR. When considering 
the correlation between metallicity and Hi mass, however, 
it can be seen that the correlation persists even when con- 
sidering the highest stellar mass bin. Even the most massive 
galaxies in our sample exhibit decreasing metallicity with 
increasing Hi mass. 

Another important point to take from Fig. [2] is the gra- 
dient of the metallicity-Hl mass slope. It can be seen that 
Hi mass correlates inversely with metallicity at all stellar 
masses, with a slope that is approximately constant with 
stellar mass. Importantly, this slope is significantly greater 
than —1 for all stellar mass bins. A slope of —1 might be 
expected if the sole physical process driving the FMR was 
ISM dilution by infalling metal-poor gas (such that increas- 
ing the Hi mass by an order of magnitude decreases the 
value of log(0/H) by a concomitant order of magnitude). 
Instead, the slope of the inverse metallicity-Hl mass correla- 
tion is closer to —0.15 (in the log-log space of Fig. 2), and is 
approximately constant across the stellar mass range probed 
by our sample. 

The fact that the metallicity dependence on Hi mass 
is essentially the same at all stellar masses suggests that 
the same mechanisms responsible for metal enrichment and 
dilution are in place across galaxies of all masses. 
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Figure 4. The edge-on projection of the Hl-FMR, which min- 
imises the scatter in metallicity. As in Fig. 1 (right), coloured lines 
show the mean trend for different Hi mass bins. The quadratic fit 
to the Hl-FMR (given in Eq. 4) has been ovcrplotted as a black 
dotted line between 8.5 < ijg < 11. 



3.2 The dispersion in the Hl-FMR 

The recent papers that defi ned the FMR (JMannucci et al] 
120101 ; lLara-Lopez et alJuOlOf ) defined the relation as a distri- 
bution of galaxies in a three dimensional parameter space de- 
fined by stellar mass, metallicity, and SFR. The three main 
'projections' of this fundamental parameter space, then, are 
the mass-metallicty relati on, the mass-SFR relat ion, and the 
SFR-metallicity relation. iMannucci et al.l (J2010j) consider a 
projection of this parameter space which minimises scatter 
and is not orthogonal to any axis (and therefore does not 
marginalise out one of the three parameters) . 

They define a variable fj, a which combines both stellar 
mass and SFR, 



ju tt = log(M), -alog(SFR), 



(2) 



and find that a value of a ~ 0.32 minimises the scatter 
in the metallicity-/^ plane. 

We take a similar approach in finding the optimal pro- 
jection of both the SFR-FMR and the Hl-FMR. For finding 
the optimal projection of the Hl-FMR, we adopt a parame- 
ter rip which combines both the stellar mass and Hi mass: 



rip = log(M), - /31og(M(Hi) - 9.80). 



(3) 



(The subtracted factor of 9.80 is equal to the (log) me- 
dian Hi mass, and simply adds a linear scaling offset to 77, 
to produce a range in 77 that is comparable to the range in 
both M» and /j.) 

We analyse the scatter in the metallicity-^ relation as 
a function of varying (3. Note that setting (3 = causes 77 
to simply equal stellar mass (so that metallicity-?^ becomes 
the standard mass-metallicity relation), and setting (3=1 
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Figure 5. Plots showing the dispersion in metallicity in the M*-M(Hl) plane (left panel) and the M*-SFRhc« plane (right panel). The 
red dotted line in the right panel shows the SFR 'Main Sequence', defined at z = by Elbaz et al. (2007). 



causes 77 to be the inverse atomic gas fraction 1//hi (if we 
define the atomic gas fraction /hi as Mhi/M*). 

Figure [3] shows the mean dispersion in the metallicity- 
H relation for SFR (left panel) and the metallicity-77 re- 
lation for Hi mass (right panel), as a function of varying 
the parameters a (left) and /3 (right). In each, the point of 
minimum dispersion signifies the optimum projection of the 
FMR. We define dispersion in two slightly different ways. We 
have taken a fourth-order polynomial best fit to the stan- 
dard a = (and j3 — 0) mass-metallicity relation following 
Mannucci et al. (2010): in each plot, the solid line shows the 
scatter around this original fit as a (and ft) are varied. As a 
check, we also re-define a fourth-order polynomial fit to the 
/i- metallicity (and 77-metallicity) relation at each separate 
value of a (and /3). The dashed line then shows the scatter 
around this fit, which varies as a function of a (and /5). It 
can be seen from Figure[3]that the two definitions of scatter 
produce minima at essentially identical points. 

In the left hand panel, we show the data for varying a, 
the proje ction of the SFR - FMR . We confirm that the result 
found by iMannucci et al.l (|2010l ) is also seen in our smaller 
sample - that the optimal projection of the SFR-FMR is 
not with a = 0, but is at a = 0.28. 

We also demonstrate, for the first time, that the FMR 
also manifests equally when considering the influence of Hi 
content on the mass-metallicity relation. The optimal pro- 
jection of the Hl-FMR is found to be with /3 = 0.35. That is, 
the metallicity of galaxies in the local Universe can be well 
described by the combination of the stellar and Hi masses, in 
the form of 770.35 = Iog(M). - 0.35 log(M(Hi) - 9.80). Fig.H 
shows the metallicity distribution of the various bins in our 
sample as a function of this new parameter 770.35, clearly 
showing that dispersion is minimised along this projection 
and that the metallicity depends solely on the parameter 
770.35- 

We have fitted the overlapping curves of metallicity ver- 
sus 770.35 with the following quadratic function: 



12 + log(0/H) = -3.979(6) + 2.363(5)770.35 - 0.1068(2)7?^ 35 , 

(4) 

which traces the Hl-FMR closely between 8.5 < 77,3 < 
11. Numbers in parentheses indicate the uncertainty in the 
final digit of the fitted coefficients. 

It is also notable from Figure[3]that the dispersion in the 
FMR varies more strongly with a (i.e. the SFR dependence) 
than of /3 (the Hi mass dependence). That is, varying the 
projection angle of the FMR has more of an effect on the 
dispersion of the SFR-FMR than the Hi-FMR. Naively, this 
may suggest that the SFR-FMR is more tightly correlated 
than the Hl-FMR (consider the extreme situation where the 
dispersion is independent of projection, and the relation in 
Figure [3] is flat; this would imply no FMR exists, and the 
parameter has no influence on the mass-metallicity relation). 
However, the optimum projection of each relation produces 
a comparable minimum dispersion: 0.073 dex for the SFR, 
and 0.075 dex for the Hi, suggesting that the FMR can be 
expressed equally in terms of both SFR and Hi content. This 
dispersion is only slightly larger than the uncertainty on the 
metallicity, which we estimate to be ~ 0.05 dex (excluding 
systematic error). 

We suggest that this slightly stronger dependence on 
SFR could be due to the non-independence of the SFR and 
metallicity diagnostics; bright emission lines (particularly 
Ha) are instrumental in calculating both the metallicity and 
SFR of the SDSS galaxies in our sample, with the result 
that these two derived parameters are not truly independent. 
(This iss ue has been noted previousl y: see lPilvugin fc Thuanl 
120051 and I Andrews fc Martimll2013l for discussion). 

Indeed, Maiolino et al. (in prep.) have found that when 
the SFR is inferred from IR emission, the dispersion around 
the SFR-FMR becomes significantly larger. This result, 
combined with the results presented in this work, could sug- 
gest that the Hl-FMR is the more fundamental relation, with 
SFR-FMR being best understood as a by-product. 
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In addition, it must be noted that this direct com- 
parison between the scatter in the SFR-FMR and the 
Hl-FMR is not entirely fair, since in the former relation 
the spectroscopically-obtained parameters (metallicity and 
SFR) are measured within the same SDSS fibre aperture, 
while in the latter relation the Hi mass is measured via in- 
tegrated 21cm emission - not from the region sampled by 
the metallicity measurements. Compounding this effect is 
the fact that HI emission is often found out to large r adii, 
beyond the stellar disk (see, e.g., ICavatte et al.lll994l ). As 
a result, it is likely that additional scatter is added to the 
Hl-FMR due to both the mismatched apertures and the ex- 
tended Hi disk. Unfortunately, we cannot confirm this di- 
rectly by extracting the Hi mass within the SDSS aperture. 
We can, however, attempt a comparison by checking the 
scatter in the SFR-FMR when using integrated SFRs (as 
inferred by Brinchmann et al. 2004), in order to introduce 
a similar 'aperture mismatch' effect to that present in the 
Hi data. In this case, the dispersion around the SFR-FMR 
is higher (~ 0.1 dex), and the relation between metallicity 
and SFR is weaker (with a closer to zero). It is likely, there- 
fore, that if we were to use matched apertures, the lowest 
scatter would be evident when considering the Hl-FMR, fur- 
ther suggesting that the Hl-FMR is the more fundamental 
relatiorQ 

Another way to visualise the scatter in these relation- 
ships (without needing to recourse to measuring the scatter 
around fitted functions) is to analyse the metallicity dis- 
persion in the M*-M(Hl) and the M*-SFR planes. This is 
shown in Figure [5] The lcr scatter in metallicity is shown 
across planes of M.-SFR (left panel) and M*-M(Hl) (right 
panel), with bin sizes as before. The inset bar shows the 
scatter colour-coding up to a maximum value of 0.3 dex. 
In order to avoid misleading dispersions at the edges of the 
distributions caused by small number statistics, only bins 
containing ^ 10 galaxies have been plotted. 

It is interesting to note that the metallicity disper- 
sion is low (around or below 0.1 dex) for intermediate stel- 
lar masses (9 < log(M»/Me) < 11) and Hi gas masses 
(log(MHi/M©) > 9.5 — 10), which may suggest that this 
region of the diagram is predominantly populated by galax- 
ies evolving smoothly, in equilibrium between SFR, gas 
inflow and outflow. In contrast, massive galaxies, with 
log(M«/Mo) ~ 11 or higher appear to have much higher 
metallicity dispersion, approaching 0.2 dex. Such a large dis- 
persion may reveal that massive galaxies are subject to a 
larger fraction of (minor) merging/interactions, which make 
the metallicity deviate from the 'regular' behaviour charac- 
terising secularly evolving galaxies. 

The region with low stellar and gas masses is also char- 
acterised by a rapidly increasing dispersion. This property of 
low mass galaxies is likely a consequence of the fact that the 
accretion of Hi clouds (responsible for the dilution of met- 
als) cannot any longer be considered as continuous smooth 
events for low mass galaxies, since the mass of individual 
in-falling clouds is likely comparable with the mass of gas 
already contained in the galaxy, hence producing a discon- 
tinuous and scattered behaviour. 



2 See Appendix A for a discussion of fibre-based aperture effects 
relevant to SDSS metallicity measurements. 



3.3 The fundamental metallicity relation and H2 
content 

Figure [5] shows the mass-metallicity relation for the 62 lo- 
cal galaxies in our sample that have been observed in CO. 
These include galaxies in the COLDGASS sample, which 
appear in the SDSS (and therefore have metallicities, stel- 
lar masses, and star formation rates calculated identically 
as for the ALFALFA sample above), as well as other galax- 
ies from the literature. These two samples are differentiated 
by plotting with full and empty symbols (respectively). In 
addition, 9 high-z SMGs with stellar masses, CO observa- 
tions, and [Nll]A6584/Ha metallicities are also plotted as 
square symbols. Both panels show the mass-metallicity re- 
lation, with different colour-coding; the left panel shows the 
mass-metallicity relation colour-coded with H2 mass, while 
the right panel colour-codes with Hi mass. 

It can be seen that there is no apparent trend in metal- 
licity between molecular gas rich and gas poor galaxies. Un- 
fortunately, due to the small sample size, separating the 
galaxies into several bins of H2 mass would not produce 
physically meaningful results. Instead, we can get a crude 
estimate of any possible trend by separating the sample into 
two broad categories: 'molecular gas rich' galaxies (with H2 
masses greater than the median value (log Mh2 = 9.35), and 
'molecular gas poor' galaxies with H2 masses lower than this. 
The mean value of the metallicity is the same for each sub- 
sample; 12 + log(0/H) = (9.00±0.09) for the H 2 rich subset, 
and (9.00 ± 0.09) for the H 2 poor subsefl The mean stellar 
masses for the two subsamples are similar, with mean stellar 
masses of log M(*) = (10.53 ± 0.38) and (10.37 ± 0.38) for 
the H2 rich and H2 poor subsets, respectively. 

This result for molecular gas content stands in contrast 
to what we have seen above (a strong dependence of metal- 
licity on Hi mass), and initially could suggest a lack of de- 
pendence of the metallicity on the molecular gas content. 
This would be a very surprising result, given that it is the 
molecular gas phase that is most intimately connected with 
star formation. However a similar lack of dependence is also 
seen in the Hi data for the H2 sample. Again adopting the 
dual 'gas rich/gas poor' binning, galaxies in our H2 sample 
which are Hi rich have mean metallicities of 9.00 ± 0.09, 
whereas galaxies in the same sample which are Hi deficient 
have a mean metallicity of 8.99 ± 0.09. The same goes for 
the relationship with the total gas masQ (9.01 ± 0.08 and 
8.98 ± 0.10 for gas rich and gas poor galaxies, respectively). 

The lack of dependence of Hi content on the metallicity 
of our smaller (H2-observed) sample, compared to the sig- 
nificant dependence observed in our large ALFALFA-SDSS 
sample (see SJ37T] above), suggests that some combination of 
sample size and selection effects are obscuring any potential 
trend. 

The COLDGASS sample (which is a random selection 
drawn from the parent sample, GASS) includes only galaxies 
with M* > 10 10 Mq, and is selected to have a flat dis- 
tribution in logM*. This stellar mass coincides with the 



3 Here, and throughout this discussion, errors represent a lcr 
standard deviation in the relevant population, rather than an un- 
certainty in the quoted mean. 

4 M(gas) = 1.36(M(Hl) + M(H 2 )), where 1.36 is a correction 
applied to account for interstellar helium. 
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Figure 6. The mass-metallicity relation for galaxies with CO observations. In each panel, galaxies appearing in COLDGASS (and 
therefore the SDSS) are plotted with filled circles. Galaxies gathered from the literature are plotted with open circles. High-2 SMGs are 
shown as squares. Left panel: galaxies colour-coded with their H2 mass. A metallicity-dcpcndent CO/H2 conversation factor has been 
assumed (see text for discussion). Right panel: galaxies colour-coded with their Hi mass. In each panel, the black point in the top left 
shows the typical uncertainty in metallicity. 



turnover in the mass-metallicity relation, with the result 
that the metallicity range of our H2 sample is restricted 
relative to our (larger) Hi sample, providing less dynamic 
range over which any potential trends can be examined. Fur- 
thermore, the 62 CO-observed galaxies remaining after our 
sample selection cuts (primarily the removal of AGNs and 
CO non-detections) exhibit a particularly restricted metal- 
licity range, with a standard deviation in 12+log(0/H) of 
< 0.1 dex. 



It is unlikely that CO non-detections are significantly 
affecting any of the results presented here. While 35% 
(83/233) of the COLDGASS parent sample are listed as CO 
non-detections, 96% of those CO non-detects (80/83) are 
removed from our sample by the selection cuts described in 
§2.1 and §2.2. Just three COLDGASS galaxies undetected 
in CO satisfy our selection criteria, and those three galaxies 
have metallicities consistent with the COLDGASS average. 
(The three CO non-detects do have a lower SFR, however, of 
log(SFR) = —1.46 ± 0.83, compared to the detected sample 
average of log(SFR) = -0.76 ± 0.62.) 



A further complication in the attempt to test the ex- 
istence of a molecular gas FMR is the strong dependence 
on the CO/H2 conversion factor on the gas-phase metal- 
licity. Although there have been several attempts to model 
and to observationally constrain the relation between the 
CO/H2 conversion factor and gas-phase metallicity (e.g. 
ICenzel et ail l201ll ; iNaravanan et all l201ll ; Bolatto et al. 
2013), such a dependence is so strong that uncertainties in 
the metallicity and/or deviations from the nominal relation 
adopted in this work, can strongly affect or wash out any 
real trend between molecular gas content, metallicity and 
stellar mass. 



4 DISCUSSION 

Several models have attempted to reproduce the Fundamen- 
tal Metall icity Relation between SFR, metallicity and stellar 
mass (e.g. lDave et al.ll2012J ; lDaval et al.ll2012l ; ICampisi et al.l 
I2011T ). in terms of inflows, outflows and star formation ef- 
ficiency. It is an important goal to investigate whether the 
same models are able to reproduce the Fundamental Metal- 
licity Relation between atomic gas mass, metallicity and 
stellar mass. An attempt to model the metalli city depen- 
dence on the Hi content was already reported in I Dave et al.l 
l|2013r ). Here we have provided a much more detailed and 
quantitative fundamental relation, which can be used to test 
galaxy evolutionary models in detail. In particular, it is im- 
portant to verify whether current models can explain the 
metallicity dependence on Hi mass even in massive galax- 
ies (in contrast to what is observed for the SFR-metallicity 
relation) and whether they can explain the large scatter ob- 
served in massive and in low mass galaxies. 

It is beyond the scope of this paper to develop detailed 
models to interpret the Hl-FMR and its scatter. We only 
note that the observed relation cannot be explained in terms 
of a simple closed-box evolutio n (as suggested by other au- 
thors, e.g. iMagrini et al.l 120121 ) where the anti-correlation 
between gas mass and metallicity is a simple consequence 
of the ageing of the galaxy (in which the gas is gradually 
consumed and enriched by star formation.) 

Fig.[7]shows the Mm vs metallicity relation for different 
stellar masses, derived using a simple 'closed box' model pre- 
sented by Erb (2008 jj. For this toy model, we have assumed 
that all gas is atomic (that is, the gas mass is equivalent 
to the Hi mass), the fraction of mass remaining locked in 



See Eq. 11 in that work. 
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stars, a = 1, and the yield y = 0.019. It must be noted that 
varying these simple assumptions does not change the form 
of the results in any significant way. We have also converted 
between 'metal content' Z (defined as the mass fraction of 
elements heavier than helium) and 12+log(0/H): 



12 + log (O/H) = log 



0.0126 



+ 8.69 



(5) 



The model data reproduce the basic form of the 
metallicity-Hl mass relation shown in Fig. 2 - that metal- 
licity is inversely correlated with Hi mass at a given stellar 
mass (and directly correlated with stellar mass at a given Hi 
mass). However, there are some significant differences which 
suggest that a simplistic 'closed box' model is insufficient 
to explain the metallicity-Hl mass relation as revealed by 
the ALFALFA sample. The model produces a metallicity- 
Hl mass correlation with a gradient that is (a) significantly 
steeper than is seen in the data, and (b) dependent on stellar 
mass. The gradient of the model metallicity-Hl correlation 
varies between —0.9 for the lowest stellar mass galaxies, to 
—0.4 for the most massive. This stands in contrast to the 
trend seen in the ALFALFA data, where the gradient is ap- 
proximately constant at —0.15 across all stellar mass bins. 
A 'pure dilution' scenario, in which metal abundance is di- 
luted by an infalling (pristine) gas cloud also fails to account 
for the data - this scenario would result in a trend between 
Hi mass and metallicity even steeper than produced by the 
closed box model. 

It is clear from these model results that the most ba- 
sic features of the data (the inverse correlation between 
Hi content and metallicity) can be reproduced by a sim- 
ple 'closed box' model, or dilution model. However, it is also 
clear that more complex modelling (likely involving both 
inflow/outflow of gas, and metal enrichment by star forma- 
tion) is required to fully account for the flatter trends re- 
vealed in Fig. 2. It is possible that outflows are responsible 
for the dependence on mass, with infalls controlling the de- 
pendence with SFR; outflows are the dominant force in the 
chemical evolution of low mass galaxies, which (in the sim- 
ple model) show the greatest discrepancy with the results 
presented here. 

In addition, further work is required to understand the 
role of molecular gas in these trends. While the CO data pre- 
sented in this work do not allow us to draw conclusions about 
any potential correlation between molecular mass, stellar 
mass, and metallicity, due to the paucity of data we cer- 
tainly have not been able to rule any such correlation out. 
Indeed, as star formation is more closely connected with 
molecular gas than atomic gas, the existence of the SFR- 
FMR would at least suggest the existence of a 'molecular 
gas FMR'. Larger datasets, and a better understanding of 
the metallicity-dependence of the CO/H2 conversion factor, 
will be required in order to explore this. 

Finally, it is an important future goal to investigate 
whether the Hi FMR discovered here also holds at high 
redshift. While hundreds of high-z galaxies have had their 
gas-phase metallicities studied, detecting Hi at cosmologi- 
cal redshifts will only be feasible with the SKA pathfinders 
(MeerKAT and ASKAP), and studying significant numbers 
of galaxies in Hi, spanning a wide range of masses, will re- 
quire the full deployment of the SKA. Potentially, the abil- 
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Figure 7. The metallicity-Hl mass relation, shown for a range 
of stellar masses, derived using the 'closed box' model given by 
Erb (2008). Compare to Fig. 2 (right panel). While the model 
reproduces some of the trends seen in the data (the inverse cor- 
relation between Hi mass and metallicity), the slope at a given 
stellar mass is far steeper in the models than is seen in the data. 



ity to investigate the relationship between metallicity and 
molecular gas will be within reach on a shorter timescale. 
Both in the local universe and at high redshift larger sam- 
ples of galaxies are expected to become available on short 
timescale thanks to ALMA, and to upcoming extensive sur- 
veys using IRAM. However, the strong dependence of the 
CO-to-H2 conversion factor on metallicity, and the associ- 
ated uncertainties, will likely remain a major issue which 
may hamper our ability to identify and characterise a puta- 
tive H 2 FMR. 



5 CONCLUSIONS 

In this work, we have examined the interconnection between 
the stellar mass, gas-phase metallicity, and gas content of 
local galaxies. Our main conclusions are as follows: 

• We have presented evidence that the mass-metallicity 
relation exhibits a strong secondary dependence on Hi 
mass, with Hi-rich galaxies being more metal poor at a 
given stellar mass. This trend is evident across the stellar 
mass range probed by our sample. 

• We introduce the Hi fundamental metallicity relation, 
an optimum projection of metallicity, stellar mass, and 
Hi mass, which reduces dispersion relative to the mass- 
metallicity relation. We expect that this is a fundamental 
relation, which lies behind the well-known SFR-FMR. 

• We also note that the dispersion in metallicity is lowest 
< 0.1 dex for galaxies of intermediate stellar and Hi mass. 
We argue that this is the regime populated by secularly 
evolving galaxies in equilibrium between SFR, inflow, and 



© 2012 RAS, MNRAS 000, 



Metallicity, gas content, and stellar mass 11 



outflow. 

• We assemble a smaller sample of galaxies observed in 
CO, but see no correlation between molecular gas content 
and metallicity. We argue that this does not rule out a 
H2-FMR, but is due to a combination of selection effects 
and complications involving a metallicity-dependent CO/H2 
conversion factor. 

• Finally, we note that a simple 'closed box' model is not 
sufficient to explain the trends in the Hl-FMR. Developing 
more sophisticated models to fully account for the behaviour 
between Hi mass, metallicity, and stellar mass is an impor- 
tant future goal. 
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APPENDIX A: APERTURE EFFECTS: 
SIZE OF THE SDSS FIBRE 



THE 



One potential source of bias that must be considered is the 
aperture effect of the SDSS fibre. Each SDSS fibre has a 
typical diameter of ~ 3". At a distance of 100 Mpc (approx- 
imately the lower limit of our Hi sample), this corresponds to 
a physical extent of a few kpc. Clearly, such small apertures 
will only capture a small fraction of the total light from more 
extended galaxies, and aperture corrections must be applied 
to photometry in order to estimate integrated magnitudes. 

The existence of abundance gradients within individ- 
ual galaxies could potentially cause metallicities calculated 
from small apertures to be biased. The galactic bulge typ- 
ically exhibits a higher metallicity than the disc (see, e.g., 
iHenrv fc Worthevlll999l ). so measuring a metallicity using 
spectra taken from a fixed central aperture can cause abun- 
dan ce estimates to be b iased high. 

iKewlev et alj 1)20051 ) estimate that a fibre covering frac- 
tion of 20% is a useful benchmark, above which abundances 
derived from a central aperture agree well with integrated 
values. We estimate the fraction of a galaxy covered by the 
SDSS fibre by taking the ratio of the fibre SFR and the total 

I i.e. ap erture-corrected) SFR, as given bv lBrinchmann et al.l 
2004h and the MPA-JHU group. Figure [5]] shows the frac- 
tion of light within the fibre plotted against stellar mass for 
our ALFALFA-Hl sample. Approximately 1/3 (1413/4253) 
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Figure Al. The fraction of each Hi galaxy contained within the 
SDSS fibre, plotted against stellar mass. The fraction within the 
fibre was calculated by taking the ratio of the fibre-SFR to the 
total, aperture-corrected SFR. Kewley et al. (2005) estimate that 
a covering fraction of > 20% (indicated with a vertical dashed 
line) is needed to avoid bias resulting from metallicity gradients. 



of the galaxies in our Hi sample have > 20% of their light 
within the SDSS fibre, and should therefore be relatively 
unaffected by bias introduced by metallicity gradients. 

We have repeated the analyses in this work for this 
'aperture unbiased' subsample. While we lose some dynamic 
range in Hi mass by restricting our sample in this way, all 
our main conclusions (the existence of a Hl-FMR, with Hl- 
rich galaxies being more metal poor at all stellar masses) 
remain robust when considering only galaxies for which bias 
due to metallicity gradients is likely to be negligible. 
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